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A muon-spin relaxation (/x + SR) investigation is presented for the molecular superconductor n- 
(BEDT-TTF)2Cu[N(CN)2Br]. Evidence is found for low-temperature phase-separation, with only 
a fraction of the sample showing a superconducting signal, even for slow cooling. Rapid cooling 
reduces the superconducting fraction still further. For the superconducting phase, the in-plane 
penetration depth is measured to be \\\ = 0.47(1) /im and evidence is seen for a vortex decoupling 
transition in applied fields above 40 mT. The magnetic fluctuations in the normal state produce 
Korringa behavior of the muon spin relaxation rate below 100 K, a precipitous drop in relaxation 
rate is seen at higher temperatures and an enhanced local spin susceptibility occurs just above T c . 
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PACS numbers: 74.70. Kn, 76.75.+i, 74.25. Ha 



I. INTRODUCTION 

The quasi two-dimensional molecular series k-(BEDT- 
TTF) 2 A exemplifies the complex interplay of collective 
phenomena that occur in correlated electron systems^ 2 .. 
In particular, the proximity of superconducting k- 
(BEDT-TTF) 2 Cu[N(CN) 2 ]Br to a Mott transition, along 
with possible pseudogap physics and molecular disorder 
effects, has led to this material being extensively studied 
in recent years^— . Despite this intense interest, several 
experimental observations in this system lack a conclu- 
sive explanation^ and more experimental study is war- 
ranted. In this paper we present a muon-spin relaxation 
(/i+SR) investigation of K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br. 
Our focus is the penetration depth in the vortex state for 
T < T c and the fluctuations of the local magnetic field in 
the normal state. 

The layered structure of the k-(BEDT-TTF) 2 X family 
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FIG. 1: (a) Layered structure of the k- (BEDT-TTF) 2 X sys- 
tem, (b) Phase diagram of the system showing the prox- 
imity of the X =Cu[N(CN) 2 ]Br salt to the Mott insula- 
tor phase. (MI=Mott Insulator, PI=Paramagnetic insulator, 
SC=Superconductor, FL=Fermi liquid.) (After Ref. @). 



of organic molecular metals is shown in Fig. [TJa). The 
BEDT-TTF molecules dimerize, forming molecular units 
which stack on a triangular lattice in two-dimensional 
planes. The removal of one electron from each BEDT- 
TTF dimer causes the tight-binding band to be half filled 
and in order to balance the charge layers of anion, X, 
are located between the partially oxidized BEDT-TTF 
sheets. The properties of the system are strongly de- 
pendent on the transfer integral t which is controlled by 
the dimer separation and may be manipulated by apply- 
ing pressure or by changing the identity of the anion 2 . 
The phase diagram of the family [Fig. QXb)] shows that 
the system may be tuned from Mott insulator through 
superconductivity into a normal metallic state as a func- 
tion of t/U, where U is the on-site Coulomb repulsion 
(which is a property of the dimer and is almost indepen- 
dent of the anion or pressure) . Importantly, this tuning 
is achieved without the need for chemical doping (in con- 
trast to the cuprates) and therefore minimizes structural 
disorder effects. At small t/U Coulomb correlations dom- 
inate and k-(BEDT-TTF) 2 A with X =Cu[N(CN) 2 ]Cl 
is a Mott insulator 8 . As pressure is increased or X is 
changed to X =Cu[N(CN) 2 ]Br there is an insulator-to- 
superconductor transition 9 . The X =Cu(NCS) 2 com- 
pound exhibits a larger t/U still and remains a super- 
conductor with a slightly depressed T c . 

Superconducting K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br 

(T c ps 12 K) displays several properties [Fig. [5] 
whose explanation remains obscure. These may be 
broadly grouped three temperature regions. In region I 
(T ps 80 K) there are sharp changes in the temperature 
dependence of the lattice constants 7 [Fig. [5Jd)] and, 
above 70 K, a rapid increase in the thermal expansion 
coefficient perpendicular to the planes^ a± [Fig. [He)]. 
These coincide with a rounded maximum in the 
resistivitjsii p [Fig. [5Jb)]. The properties of this region 
have been linked to a glass-like freezing of terminal 
ethylene groups on the BEDT-TTF molecules^ around 
T g = 77 K, although is has been suggested^ that this 
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FIG. 2: Comparison of several experimental properties with 
regions I (T » 80 K), II (T w 50 A") and III (T c < T < 
18 K) shaded, (a) LF /i + SR relaxation rates (this work); (b) 
Resistivity^; (c) NMR 1/TiT— ; (d) lattice constants^; (e) 
thermal expansion coefficient perpendicular to the layers^. 



This paper is structured as follows: after discussing 
experimental details in section II we present the results 
of transverse field muon-spin relaxation measurements 
of the superconducting penetration depth of k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br in section III. Section IV describes 
the use of longitudinal field measurements to investigate 
spin fluctuations in the normal state of the material and 
our conclusions are presented in section V. 



II. EXPERIMENTAL 

Muon-spin spectroscopy-^ is a sensitive means of 
probing the superconductivity and local magnetism of 
molecular materials of this sor t 14 i 15 . Transverse-field 
(TF) and longitudinal-field (LF) muon-spin relaxation 
(/i + SR) measurements were made on a mosaic sam- 
ple of K-(BEDT-TTF) 2 Cu[N(CN) 2 Br] at the ISIS facil- 
ity, Rutherford Appleton Laboratory, UK and the Swiss 
Muon Source (S/xS), Paul Scherrer Institut, CH. The 
polycrystalline sample was made of ~ 50 small crystal- 
lites which grow as platelets whose large faces are parallel 
to the conducting layers (the ac planes for this material). 
These were arranged on a piece of Ag foil to cover an area 
of w 0.5 cm 2 so that the conducting layers are parallel 
to the plane of the mosaic. TF measurements were made 
using the MuSR spectrometer at ISIS where the sample 
was mounted on a hematite backing plate in a helium 
cryostat with the sample oriented at 45° to both the ap- 
plied magnetic field and the initial muon spin. In order 
to reduce any effect of ethylene disorder on the supercon- 
ducting properties, the sample was slow cooled at a rate 
of < 5 K/hour. TF measurements were also made using 
the GPS instrument at S/iS where the sample plane was 
oriented at 90° to the applied magnetic field and no slow 
cooling procedure was followed. LF measurements were 
made on the ARGUS spectrometer (ISIS) in the latter 
90° geometry. 



III. TF ^+SR AND SUPERCONDUCTIVITY 



freezing alone cannot account for the observed structural 
changes. In region II (T ss 50 K) NMR measurements^ 1 ^ 
yield a maximum in 1/TyT around T* w 50 K [Fig. fS^c)] . 
Upon cooling through this region there is a crossover 
in transport properties from "bad metal" behavior to a 
more conventional Fermi-liquid regime 2 - and a plateau in 
a± is also seen [Fig. f2Je)]. These phenomena have been 
linked to pseudogap physics although the details remain 
unclear 2 . In region III (T c < T < 18 K) a vortex Nernst 
signal is observed^. It was argued- that the proximity of 
the X =Cu[N(CN) 2 ]Br material to the Mott state in the 
phase diagram in Fig.fTJb), results in vortex fluctuations 
persisting above T c and, furthermore, this remnant of 
fluctuating superconductivity may be consistent with the 
occurrence of phase fluctuations in the superconducting 
order parameter close to the Mott boundary. 



TF /^ + SR provides a means of accurately measuring 
the internal magnetic field distribution in a material, 
such as that due to the vortex lattice (VL) in a type 
II superconductor 1 -. In a TF /i + SR experiment spin po- 
larized muons are implanted in the bulk of a material in 
the presence of a magnetic field B c \ < B a < B c2 , which 
is applied perpendicular to the initial muon spin direc- 
tion. Muons stop at random positions on the length scale 
of the VL where they precess about the total local mag- 
netic field B at the muon site with frequency w^ = 7 M -B, 
where 7^ = 2ir x 135.5 MHz T _1 . The observed property 
of the experiment is the time evolution of the muon-spin 
polarization P x {t), which allows the determination of the 
distribution p(B) of local magnetic fields across the sam- 
ple volume via P x (t) = J °° p(B) cos^^Bt + 4>)dB, where 
the phase <f> results from the detector geometry. 
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FIG. 3: (a) Broadening of the 1.5 K TF /i+SR signal due to 
the VL measured at S/iS. The field B a was applied perpen- 
dicular to the conducting layers. The dotted line indicates 
a transition to a decoupled-layer vortex phase (shaded re- 
gion) occurring between 40 mT and 80 mT. The solid line 
is the fit of the field dependent line width 1 ^ below the vor- 
tex transition, giving Xu = 0.47(1) /jm; the dashed line is the 
corresponding width predicted by Eq. {TJ. The dot-dash line 
shows the width according to Eq. ([2]) for fully decoupled vor- 
tex layers, (b) Temperature dependence of the normalised VL 
broadening measured at ISIS in 3 mT field applied at 45° to 
the sample plane. A fit is shown to Eq. J3J. The inset shows 
scaling 2 -^ between T c and A7 2 , comparing the current data 

(Br) with that of k-(BEDT-TTF) 2 Cu(NCS) 2 (NCS)^; the 
open circle shifted to the left shows the apparent effect of 
applying B a > 80 mT to Br. 



Above T c some broadening of the spectrum is caused 
by randomly directed nuclear moments relaxing the 
muon spins. Below T c the spectrum broadens consider- 
ably due to the contribution from the VL, but this com- 
ponent was found to be only ~ 30% of the total signal 
for T < T c . Allowing for part of the non-superconducting 
signal to come from the cryostat and sample mounting 
material, we estimate that at most 50% of the sam- 
ple is giving a superconducting signal. This is signif- 
icant since it suggests that the majority of muons are 
not stopped in regions where there is a well defined VL 
and may indicate the coexistence of a competing phase 
in this material. In fact, there is evidence for a phase 
separation between superconducting and insulating re- 
gions in this material^ which is strongly dependent on 



disorder, such as that introduced through rapid cooling. 
The non-superconducting contribution to our signal for 
T < T c decreases weakly and linearly with increasing 
temperature. Rapid cooling of the sample was found to 
reduce the superconducting volume fraction still further. 
With such a small superconducting fraction it is diffi- 
cult to accurately extract detailed lineshape information 
from our data, however the RMS width of the VL field 
distribution -Bvl can be straightforwardly obtained by 
fitting the data to the sum of two Gaussian relaxation 
components, reflecting the superconducting VL and non- 
superconducting contributions. These are easily sepa- 
rated below T c since the contribution from the VL has 
a much larger relaxation rate. The field dependence of 
the VL component is shown in Fig. I^a). The low field 
data (40 mT and below) fit well to the width expected 
from the solution of a Ginzburg-Landau (GL) model of 
the triangular VL field distribution— with an in-plane 
penetration depth of 0.47(1) /im. -Bvl is not very sensi- 
tive to B C 2 in this field range, so we assume the reported 
value 18 B C 2=10(2) T; the GL parameter n is then esti- 
mated to be 80(10). The penetration depth in the highly 
conducting planes can also be estimated via the approx- 
imate expression for a high anisotropy superconductor 
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where $o is the flux quantum and 9 is the angle of the 
applied field with respect to the normal to the plane. 
This field-independent expression is valid at intermediate 
fields B c i<^B a <^B c2 and the -Bvl value obtained from 
Eq. ([I]) for An = 0.47 /im is shown as the dashed line 
in Fig. Efa); this is seen to be reasonably close to the 
field-dependent GL width in a field region centred around 
10 mT. 

Our value of Am is rather smaller than previous re- 
ports, e.g. the range 0.57 < A|| < 0.69 /im obtained from 
reversible high field magnetization measurements^, and 
it is significantly lower than the value Aii ~ 0.78 /im es- 
timated from a previous /i + SR measurement^ in a field 
of 300 mT. These differences can be understood by not- 
ing that the measured linewidth decreases sharply for 
B & > 80 mT (Fig. [^a)). If this suppressed linewidth 
were taken to represent a full 3D VL, a penetration depth 
of order A = 0.80 /im would also be obtained from our 
data using Eq. ([T]) , in good agreement with the previous 
/i + SR result. The decrease in Bvl with applied field that 
we observe is not surprising, given the vortex phase dia- 
gram in (BEDT-TTF) 2 Cu(SCN)2^2i. In that material 
the ideal 3D triangular VL which exists at low tempera- 
tures and low fields^ is destroyed by the application of 
fields above a disorder-dependent threshold, that can be 
as low as 6 mT, causing a transition to a decoupled-layer 
vortex glass phase, accompanied by a sharp decrease in 
the measured linewidth—. It is therefore likely that a 
similar transition to a decoupled-layer vortex glass phase 
occurs in the X =Cu[N(CN)2Br material above a thesh- 
old field in the 40-80 mT region. The effect on the width 



TABLE I: Comparison of penetration depth Ay extracted from /i + SR measurements under different experimental conditions. 
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of the field distribution in the limit of losing all inter-layer 
correlation has been calculate d 23 ' 24 , giving 



B 2 u/B YL = lAis/a) 1 ' 2 , 



(2) 



where s is the inter-layer spacing and a is vortex spacing 
within the layers. This limit is shown as the dot-dash line 
in Fig. EJa) and it can be concluded that some residual 
interlayer correlations remain here since the data points 
lie significantly above this limit. A further feature of our 
data in the decoupled vortex phase is that the supercon- 
ducting fraction recovers to the larger value obtained in 
measuremants made under slow-cooling conditions, sug- 
gesting that the disorder-induced phase separation is be- 
coming significantly reduced once the layers are decou- 
pled. 

The temperature dependence of -Bvl has been mea- 
sured in a field of 3 mT applied in a 45° geometry and 
this is shown in Fig. [3Jb) . This can be fitted to an em- 
pirical power law 



s VL (T) = B VL (o)[i-(r/r c n, 



(3) 



where -Bvl(O) is the zero temperature contribution from 
the VL. The VL here therefore appears to thermally sta- 
ble, in contrast to the case of (BEDT-TTF) 2 Cu(SCN) 2 
where a clear melting of the VL is seen 2 ^. The val- 
ues B V l(0)=0.31(1) mT and T c =11.3(4) K are obtained 
from the fit. The penetration depth in the highly con- 
ducting planes could be obtained by taking 9 — 45° in 
Eq. {1}, however the absolute accuracy is reduced com- 
pared to 90° measurements for two reasons: firstly the 
width measured at 45° is particularly sensitive to any 
angular misalignment between the sample and field and 
secondly the angular scaling in highly anisotropic super- 
conductors can be much more complex than the simple 
form implied by Eq. $T$£. Returning to the 90° data, 
we note that a marginally better fit is obtained by allow- 
ing for an asymmetric lineshape 14 parameterized using 
the skewness parameter /3 = ((B) — B a )/BvL- The op- 
timum fit is achieved with f3 = 0.7, which is reasonably 
close to the ideal triangular lattice value of 0.60, that was 
observed in more detailed studies of the vortex phases 
of the X=Cu(NCS) 2 compound 21 -. This fit also yields 
0.47(1) fim for A||. The values obtained here and in pre- 
vious muon studies for 90° geometry are summarized in 
Table HI We may compare our values of Aii and T c for 
the X=Cu[N(CN) 2 ]Br material to values obtained previ- 
ously for the X=Cu(SCN) 2 material. This is shown as 



the inset to Fig. [3Jb) and suggests that there is a scaling 
relation between these two parameters. Although this 
pair of K-phase points lie rather close, a simple linear 
scaling 27 would be reasonably consistent with the data. 
More generally, a T c oc A~ 3 scaling has been suggested to 
apply when a broader range of molecular superconduc- 
tors is considered 2 ^, but more data would be required 
to test whether this holds or not within the subset of 
K-phase BEDT-TTF compounds. 



IV. LF /i+SR AND SPIN FLUCTUATIONS 

We now turn to the local magnetic fluctuations in the 
normal state of this system, which we probe with LF 
/i + SR measurements made in a field of B = 2 mT. In 
these measurements the initial muon spin is directed par- 
allel to the applied field and we measure the polarization 
P z (t) along the same direction via the muon asymmetry 
function A(t)[oc P z (t)]. Dynamics in the local magnetic 
field distribution will cause muon spins to flip leading to a 
depolarization with relaxation rate A. The small applied 
longitudinal field was intended to quench the contribu- 
tion to the spectra of background from nuclear moments 
allowing the contribution due to electron spin dynamics 
to be discerned. 

A typical LF /i + SR spectrum measured at 20 K is 
shown in Fig. (Ha). The relaxing amplitude is quite small, 
but this is not unexpected given previous zero field muon 
results measured on a related BEDT-TTF material^. In 
that study, muons were shown to be sensitive to mag- 
netic order within the organic layers but only a small 
relaxing asymmetry was observed. We therefore expect 
that our results will reflect magnetic fluctuations within 
the organic layers of k-(BEDT-TTF) 2 V. Above 100 K 
our spectra are best described with a single exponential 
function with relaxation rate A characteristic of dynamic 
fluctuations in the local magnetic field at the muon sites 
in the material. For T < 100 K, there is a single, highly 
damped, temperature dependent oscillation visible for 
times t < 4 fj,s. In this context, spectra of this form 
suggest a contribution from quasistatic disordered mag- 
netic moments which is typically described by a Kubo- 
Toyabe (KT) function^ / KT (A, B, t) where the data are 
best fitted with a constant field width A = 0.72 MHz for 
all T < 100 K. We note that this is quite large for ran- 
domized nuclear fields (where typically A < 0.3 MHz) 
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FIG. 4: (a) Typical spectra measured in a LF of B — 2 mT at 
T = 20 K. (b) Temperature evolution of A/T. Regions I, II 
and III are shaded and a fit to a model of antiferromagnetic 
spin fluctuations is shown. Inset: Evolution of A with T 
extracted from fits to Eq. Q. 



and given the small size of the electronic moments ob- 
served with muons in BEDT-TTF layers previously^ it 
may be that this signal results from a population of static 
disordered electronic moments. The origin of such qua- 
sistatic electronic moments is unclear although we note 
that there is no decrease in A with increasing T as would 
occur in the presence of magnetic order. We note that it 
is possible that this signal arises from those portions of 
the sample not contributing to the superconducting VL 
in TF measurements. 

The dynamic contribution to the spectra for T < 100 K 
was found to be best modelled by multiplying the KT 
function by an exponential factor e * and the data were 
fitted to the resulting functional form 



A{t)=AofK T (A,B,t)e 
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A h 



(4) 



The extracted relaxation rates A are shown in Fig. [Ha) 
and inset in Fig. HJb). Other parameterizations are pos- 
sible but lead to less successful fits. However, for com- 
parison, a fit across the measured temperature regime of 
a single exponential with relaxation rate A is shown in 
Fig. [H[a), where we see that both track the same behav- 
ior. 

At temperatures T < 100 K, the dominant trend is a 
roughly linear increase in relaxation rate with T up to 
around 100 K, above which it drops precipitously [inset 
Fig. HJb)]. The linear regime observed for T < 100 K 
is typical of a Korringa law relaxation^, resulting from 
flip-flop transitions of electronic and muon spins, and 
which is expected to lead to a relaxation rate A = ^- oc 



k BTY^ q A(q) 2 lim UJ ^ox"(q,u)/uj, where x"(g,w) is the 
imaginary part of the dynamic magnetic susceptibility 
and A(q) is the hyperfine coupling. For the case of a 
q-independent coupling, we obtain A oc fcBT[Ax(0,0)] 2 , 
showing that the muon relaxation rate probes the local 
spin susceptibility x(Q — 0> w = 0). The sharp decrease 
in relaxation rate around 100 K broadly coincides with re- 
gion I, which may suggest a link with the ethylene disor- 
der effects and glass-like transition in this region. Given 
the Korringa law behavior observed at low temperatures, 
it is likely that the crossover to nonmetallic behavior of 
the system near region I [Fig. HJb)] changes the nature 
of the muons' coupling to the material and that at these 
elevated temperatures a different, dominant relaxation 
channel presumably opens up. To further investigate the 
temperature evolution of the local susceptibility from the 
Korringa relaxation we plot A/T(= 1/TT) against T in 
Fig. [4fb). We see that on warming above T c the quan- 
tity A/T peaks at around 15 K before decreasing sharply 
around region III (15 < T < 25 K) indicating a sig- 
nificant decrease of the local magnetic susceptibility in 
this region. Above 30 K A/T decreases only gradually 
throughout region II until dropping off suddenly in re- 
gion I. 

Our muon A/T results are quite different to the NMR 
1/TT results [Figf5] where a parger peak is seen in re- 
gion II. We do not see any dramatic change of behavior 
in region II in the /x + SR results, while there is no sug- 
gestion of any anomaly in the NMR behavior in regions I 
and III. In fact, the NMR results for 50 K< T < 300 K 
are well described by a model of antiferromagnetic spin 
fluctuations^ based on the MMP model 3 - , which predicts 
1/TT ra A + B/(T/T x + l), where A and B are constants 
that depend on the correlation length of the spin fluctu- 
ations whose energy scale is determined by the tempera- 
ture parameter T x . As shown Fig.^Jb) the model fits the 
muon A only passably well for T < 100 K becomes gradu- 
ally worse as region III approached from above. However, 
given that the muon sites in this material are likely to be 
different to the 13 C site where the nuclear resonance is 
achieved^, it is possible that the muon will probe a quite 
different field distribution to that probed in NMR. We 
note further that our muon results are also quite differ- 
ent to the LF /i + SR behavior observed in molecular su- 
perconductors such as the alkali-fullerides 33 where A/T 
is quite featureless above T c . 



V. CONCLUSIONS 

In conclusion, our investigation of the supercon- 
ducting and normal state properties of k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br using implanted muons has allowed 
us to identify a number of experimental features of this 
system. In the superconducting phase two important ef- 
fects are seen: firstly, a reduced superconducting signal 
fraction that is affected by sample cooling rate and is 
consistent with the phase separation suggested by ear- 



lier studies and secondly, a vortex transition taking place 
between 40 mT and 80 mT, that reduces the width of 
the magnetic field distribution for higher measurement 
fields. By taking these two effects carefully into account 
we have obtained an improved and more reliable estimate 
for the T = in-plane penetration depth Aii =0.47(1) /iin 
and find that the trend of T c increasing with superfluid 
stiffness p s ex 1/Ajj for k-(BEDT-TTF) 2 X superconduc- 
tors is consistent with the overall trend for molecular 
superconductors, whereas previous data had suggested 
that K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br was anomalous in 
this regard. In the normal state we find a large peak in 
the longitudinal muon spin relaxation rate around 100 K, 
coinciding with the region where ethylene disorder effects 
are prevalent. We also observe an enhancement of the 
local magnetic spin susceptibility above T c where vortex 



fluctuations occur in Nernst measurements. It would be 
desirable in future to extend /^ + SR studies of the normal 
state magnetic properties to the other members of this 
series in order to follow how these features evolve across 
the phase diagram. 
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